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ABSTRACT 

Recently iron emission features have been observed in several X-ray afterglows 
of GRBs. It is found that the energy obtained from the illuminating continuum which 
produces the emission lines is much higher than that of the main burst. The observation 
of SN-GRB association indicates a fallback disk should be formed after the supernovae 
explosion. The disk is optically thick and advection-dominated and dense. We suggest 
that the delayed injection energy after the initial main burst, much higher than energy 
of the main burst, causes the re-burst appearance in GRB afterglow and illuminates 
the region of the disk surface with r « l(r is the optical depth for the Thomson scatter 
) and produces the iron emission line whose luminosity can be up to 10 45 ergs" 1 . The 
duration of the iron line emission can be 10 4 - 10 5 s. This model can explain the 
appearance of re-burst and emission lines in GRB afterglow and disappearance of the 
iron emission lines, and also can naturally solve the problem of higher energy of the 
illuminating continuum than that of the main burst. This scenario is different from the 
models put forward to explain the emission lines before, that can be tested by SWIFT 
satellite. 

Subject headings: gamma rays:bursts-line:profiles-accretions disks-supernovae:general 



1. Introduction 

Recently evidence is mounting that 'long duration' GRBs are associated with a rare type of su- 
pernovae event such as a 'failed supernova' ,'Hypernova' or 'Collapsar'( Woosley 1993; Paczyhski 
1998; MacFadyen & Woosley 1999; MacFadyen,Woosley and Herger 2001; Proga et al. 2003). In 
this model, the time between the supernovae explosion and the GRB is very short, nearly simulta- 
neous. Another model called 'Supranova model' was proposed by Vietri & Stella (1998), in which 
the supernova(SN) explosion initially results in the formation of a comparatively massive, mag- 
netized neutron star endowed with rapid rotation. This supramassive neutron star is envisioned to 
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gradually lose rotational support through a pulsar-type wind until it eventually becomes unstable to 
gravitational collapse,leading to the formation of a BH and the triggering of a GRB. In 'Supranova 
model', the time between the supernovae explosion and the GRB can range from several weeks to 
several years. 

After the supernovae explosion, the fallback material will form a disk around the newly central 
black hole in either 'collapsar model' or 'supranova model', and this disk is advection-dominant 
and hot and dense even after the GRBs(e.g. Chevalier 1989; Mineshige et al. 1997). 

X-ray emission lines observed in X-ray afterglow of GRBs provide important clues for iden- 
tifying the nature of the progenitors of long (t>2 s) GRBs. The first marginal detection of an 
emission line was in the X-ray afterglow of GRB 970508 with the BeppoSAX NFI (Piro et al. 
1999). Later emission lines were also detected in the X-ray afterglows of GRB 970828 (Yoshida 
et al. 2001) with ASCA; GRB 991216 (Piro et al. 2000) and GRB020813 (Butler et al. 2003) with 
Chandra; GRB 011211 (Reeves et al. 2002), GRB 00 1025 A (Watson et al. 2002) and GRB 030227 
(Watson et al. 2003) with XMM-Newton; GRB 000214 (Antonelli et al. 2000) with BeppoSAX. 
The detailed properties of the X-ray emission features can be found in several papers ( Lazzati 
2002; Bottcher 2003; Gao & Wei 2004 ). The energy of the emission line found in the X-ray 
afterglows of GRB 970508, GRB 970828, GRB 991216 and GRB 000214 is roughly consistent 
with Fe K a at the redshift of the hosts. It is adduced as evidence that the environment of the burst 
is heavily enriched in iron as the result of a recent supernova explosions ( e.g. Lazzati et al. 1999; 
Ghisellini et al. 1999). We would have observed Fe line if the time delay is more than several 
months such as models of Supranova (Vietri & Stella 1998). 

The standard model for gamma-ray burst afterglows assumes that relativistic material is de- 
celerating on account of interaction with the surrounding medium,with a nearly impulsive injection 
energy. But perhaps the ejecta consists of many concentric shells moving at different speeds, slow 
moving material carries most of the system's energy. The delayed injection energy can be more 
than that of the main burst. It is proposed to be refreshed shock scenario. In the refreshed shock 
scenario, assuming the source ejects a range of Lorentz factors with the mass M(> 7) oc 7" v (Rees 
& Meszaros 1998; Sari & Meszaros 2000). 

The re-burst phenomena have been found in the afterglow of GRB970508(Piro et al. 1998) 
and GRB970828(Yoshida et al. 1999). It has been considered that energy of delayed injection with 
more than that of the initial burst caused the re-burst appearance(Panaitescu et al. 1998; Kumar & 
Piran 2000; Sari & Meszoros 2000). 

The energy contained in the illuminating continuum that is responsible for the line production 
is much higher than that of the collimated GRBs (Lazzati 2002; Ghisellini et al. 2002; Gao & 
Wei 2004). In this letter we consider the illuminating continuum which produces the emission line 
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comes from the delayed injection energy , which illuminates the fallback disk and produces the 
emission lines. 



2. Delayed energy injection and Fe emission line feature in X-ray afterglow 

The re-burst of emission during the afterglow has been reported in two GRBs,GRB970508 
(Piro et al. 1998) and GRB970828 (Yoshida et al. 1999). At the same time possible evidence for 
the existence of Fe Ka line has also been found in these two GRBs (GRB970508, Piro et al.1999; 
GRB970828, Yoshida et al. 1999). Delayed energy injection (or refreshed shock model) can well 
explained the resurgence emission in these two GRBs (Panaitescu et al. 1998; Kumar & Piran 
2000). The delayed injection energy more than the initial fireball energy can produce the observed 
re-burst in the afterglow about 0.5 day after 7-ray burst. Here we model that the delayed injection 
energy illuminates the fallback disk around the central black hole and photonionizes the layer of 
the disk with r=l, and iron line emission can be produced by the recombination process. 

Considering an engine that emits both an initial impulsive energy input as well as a contin- 
uous luminosity, the latter varying as a power law in the emission time. The differential energy 
conservation relation in the observer's frame can be expressed as (Cohen & Piran 1999; Zhang & 
Meszaros 2001) 

dE/dt® = L^/hf-kiE/t^.t® > t . (1) 

The first term on the right side is L=L (%/?o)^ which is the intrinsic luminosity of refresh 
shock, to is a characteristic timescale for the formation of a self-similar solution, E and t® denote 
the energy and time measured in the observer frame, q and k are dimensionless constants. 

In the refreshed shock scenario, following Rees & Me'szaros (1998), Kumar & Piran(2000), 
Sari & Meszaros (2000), one can obtain the relationship between the temporal index a and the 
spectral index (3, where F u oct^is 13 . For the X-ray afterglow, considering the forward shock in the 
slow-cooling regime(Sari, Piran, & Narayan 1998), a=(l-q/2)/3+l+q (Zhang & Meszaros,2001). 
Since no spectral information was available for GRB970508 in the first several days of the after- 
glow, in the calculation of Panaitescu et al. (1998), s=1.5 is needed. For the forward shock in the 
slow-cooling regime in the refreshed model, a = - 6 ~ 6 ( ^ /3 (Sari & Meszaros 2000, for the uniform 
medium environment). So /?=(17a-3)/24. In the X-ray afterglow observations of GRB970508, the 
temporal index a changes from -1.1 (before the re-burst) to +1.7 (at the beginning of the re-burst) 
to -0.4 (during the re-bust) and then to -2.2 (after the re-burst) (Piro et al. 1999). 

From above in the refreshed shock scenario one can obtain the delayed energy injection as 
Loc t®, q~-0.8, which ends after about 10 5 s. We assume the delayed energy is isotropic, that is 
reasonable from the observations. On the observational side, Pedersen et al.(1998) have found 



-4- 



the optical light curve of GRB970508 afterglow can be explained in terms of an isotropic outflow. 
The radio observation of GRB991216 afterglow also showed there was an isotropically energetic 
firebail(10 54 erg)(Frail et al. 2000). From the numerical work of Panaitescu et al.(1998), the injec- 
tion energy injects into the GRB outflow up to be E in j = 3E into the external medium, whose solid 
angle is fly, f2 7 is the solid angle of the GRB collimated jet. We name this energy as Injection 
Energy, that is only a fraction of the whole delayed energy. The whole delayed energy is about 
Edei = i^-E in j. Luminosity varies as about ?~ - 8 . At the beginning of the re-burst, t t ~ 6 x 10 4 s, the 
Injection Energy has been as much as that of the initial burst £ ~4x 10 50 ergs(Bloom et al. 2003). 
After that time, the residual Injection Energy 2E Q has been exhausted in about 10 5 s(Piro et al. 
1999). From 



Ju 



Li(-)-°- 8 dt (2) 



t t means the time at which the re-burst appears, t e means the time when the re-burst ends, L, 
is the luminosity at the time of re-burst emergence, E means the energy injecting into the external 
medium. After the emergence of the re-burst, the residual Injection Energy is about 2Eq, and time 
duration is about 10 5 s, we can obtain that L, ~ 10 46 ergs' 1 . From L oc r° 8 and about energy of E 
has been exhausted before the re-burst, we can get L ~ 10 41 ergs' 1 , at t ~ 3 x 10 3 s. 

An advection-dominant disk can exist around a stellar black hole even after the 7-ray bursts(Kohri 
& Mineshige 2002;Janiuk et al 2004)(Fig.l). For the 'collapsar model', SN explodes at almost 
same time with the GRBs, or about a minute to few hours prior to the 7-ray bursts. Whereas for 
'supranova model', the time delay between the SN explosion and GRB is perhaps several months 
or even longer. In this case almost all the fallback nickel has been decayed to iron. 

The evolution of fallback disk around the black hole has been considered by, eg., Meyer 
& Meyer-Hofmeister(1989), Cannizzo et al.(1990), Mineshige et al.(1993), and Mineshige et 
al.(1997). From the work of Mineshige et al. (1997),we can draw the accretion rate of the disk 
result from the fallback material. Since the radioactivity time scale is about 85 days for nickel 
decaying to iron, here we adopt 'Supranova model'(Vietri & Stella 1998), and assume 7-ray bursts 
take place about 100 days after the SN explosion, . 

M ~ 10^-'(^)(-^)(-^)(-i-r"= (3) 
3M Q 0.1M© 0.01 lOOday 

where AM is the amount of fallback material of the disk and a is the viscosity parameter. We adopt 
that the mass of the black hole to be 3 Af . AM has a large range for the different mechanism of 
the SN explosion and the evolution of the disk (eg. Woosley 1993; Chevalier 1989). Here we adopt 
AM to be 0.1 M Q . 

While about one hundred days after the SN explosion, the accretion rate would be decreasing 
to about 1O" 8 M 5 _1 . At this time the radiation pressure is dominant in the disk (eg. Mineshige et 
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al.1997). The relations of T-E and M-E are as following (Kohri & Mineshige 2002): 



T = 4.87 x 10 1 a°Eir _ 5Af| H (4) 



, i i 



~\3»>^-l„-i 



E = 3.3 x WMa'r-m B 2 H (5) 
The total mass of the disk is 

AM = / lirRHdR (6) 

J Rin 

Ri n and R represent the innermost radius and the outermost radius respectively. From eq.(3), 
eq.(5) and eq.(6), adopting AM=O.1M , we can obtain the radius of the outermost disk, R ~ 
2 x 10 12 cm. 

So we can get the results about surface density E and temperature T at the outermost radius: 

E = 4. 5 xioV».- 2 ( T ^ 7 r)(^r' (7) 

, r i M BH i 

X( 7—) 2 ( ) 2 



and 

r = 2x 10 6 iS:(-— p(— — )? (8) 
v 10 6 r, V 3M0 v 

Then the average density of the disk at this time at the outermost radius r=10 6 r v is: 

p = ^ = 7.2xl0-V^( 45x ^ Vm _ 2 ) 2 (9) 

x( T )-4 ( r 3 M 
V 2xl0 6 iT V 10 6 r/ V 3M ; 

At the outermost radius which is about r=10 6 r s =8. 85 xlO n cm, temperature T is about2xl0 6 K 
and the number density n is about 4x 10 19 cm" 3 in the disk. 

The ionization parameter is f = = 2.5 x 10 3 ( w ^- r )(4^V)" 1 ( T o^)" 2 - At this ioniza- 
tion parameter, iron emission is very efficient(Lazzati et al. 2002). The recombination time for 
hydrogenic iron in the outer disk photoionized by the nonthermal delayed energy is(Lazzati et al. 
1999) 

t rec = 1.5 x lO^rjn^s (10) 
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The temperature parametrization used here is consisitent with the range expected from photoion- 
ization equilibrium (Lazzati et al. 1999). 

The optical depth at outer radius r=10 12 cm is optically thick: T T =2nHa T ~ 10 7 . So the number 
of Fe nuclei in the layer of the disk with r = 1 is N Fe ~XFeM./(T T 56m p )~10 47 XF e - (XFe is the iron 
mass fraction of the disk ). The Fe line luminosity is [N fe (8 KeV)/t rec ](l+z) _1 , or 

L Fe ~ XFe \tf\\ +z y'ergs- 1 (11) 

For SN1987A, XFe can be about 2% when all the nickel have decayed to iron(Chevalier 1989). 
So the luminosity of iron line can be obtained: L Fe ~2x 10 45 (l+z) _1 erg5 _1 . 

After the emergence of the re-burst, luminosity decays from about I0 46 ergs~ l at the rate of 
t~° s . So the luminosity of the iron line should decrease and disappear during the re-burst, consistent 
with the observations of iron line(Piro et al. 1999). 

It is noticed that we assume the delayed energy is almost isotropic. The energy obtained by 
the disk is E disk = ^E deh Vl d is the solid angle subtended by the fallback disk as observed at the 
location of central engine. For the advection disk, H/R 0.77(Narayan et al. 2001). We can get 
E d « 031E del . For GRB970508, the open half angle of the GRB collimatedjet is 16.7°(Frail et al. 
2001). So the energy obtained by the disk is E d ~ 60E . Even if only 10% fraction of the energy 
was reflected by the disk(e.g. Zycki et al. 1994), It is about 3E reflected by one surface of the 
disk, that is sufficient for the line emission production. 



3. Discussion and Conclusions 

It is found that the energy contained in the illuminating continuum which is responsible for 
the line production is much higher than that of the collimated main GRBs. Here we model that the 
energy obtained from the delayed injection energy, higher than that of the collimated 7-ray bursts, 
illuminates the fallback disk which is formed after the supernovae explosion, photoionizes the disk 
region of r = 1, then produces the observed iron line feature. 

In our model the delayed energy comes from the central engine, which can be the mag- 
netic energy from the declining magnetic field of the superpulsar(Rees & Meszaros 2000) 
or the magnetic dipole radiation of the magnetar(Dai & Lu 1998). It could be primarily in 
a magnetically driven relativistic wind(which would be super-Eddington). The magnetized 
wind would develop a shock before encountering the disk. The nonthermal electrons would 
be accelerated behind the shock in the outflow material. The shock-accelerated electrons 
could cool promptly, and would yield a power-law X-ray continuum. It is similar to what 
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has been proposed by Rees and Mesz<iros(2000). This X-ray continuum illuminates the fall- 
back disk and produces the iron line. The surface of the disk can be accelerated outward by 
this super-Eddington flux of the illuminating continuum(e.g. Vietri et al. 2001), so usually a 
outward velocity can be seen in the lines(e.g. Reeves et al. 2002). In our model, the delayed 
energy emission and the GRBs could come from different physical processes. The delayed en- 
ergy could be from magnetic wind of the magnetar, so the energy emission is almost isotropic. 
However, from the observation of the GRBs, the GRB prompt emission should be intrinsi- 
cally collimated. So there should exist a transition from a collimated to an uncollimated 
energy release in the engine. 

The Injection Energy must be higher than that of the initial main burst so that the effect can be 
observed in the GRB afterglows(Cohen & Piran 1999;Zhang & Meszaros 2001). For GRB970508, 
at t ~ 3 x 10 3 s after the GRB, the delayed illuminating continuum decays to be 10 47 ergs' 1 , and 
the ionized iron ion recombines, Fe line appears and the line luminosity is about l0 44 -\0 45 ergs' 1 . 
at the time t ~ 6 x 10 4 s, the re-burst emerges. And after that time, the luminosity of the delayed 
illuminating continuum decays to be less than \0 46 ergs' 1 , the iron line would decrease and dis- 
appear. The duration of the Fe line tj ~ 10 4 - 10 5 s, longer than the cooling time of thermal disk 
tcooi ~ 10" 4 n7gS. All above are consistent with the observations of the Fe line in GRB970508 X-ray 
afterglow. 

The re-burst phenomenon has not been observed in GRB99 1216 afterglow. It can be explained 
that the Injection Energy is as much as or less than that of the main burst, or the re-burst was 
missed in the observation though it had happened. In the former case, the line duration should 
be 10 4 s or less, consistent with what has been observed, adopting the energy of the main burst 
E~ 10 51 ergs(Bloom et al. 2003). 

Above scenario is based on the supranova model in which the time delay between the SN 
explosion and the 7-ray bursts can be several months or even longer. Our model supports 
the 'supranova model' because it must have enough time to let nickel decay to iron (about 
10~ 3 M o Fe in the disk). In our model we assume the time between the SN explosion and the 
GRB is about 100 days, so the fallback disk that we consider has evolved for about 100 days 
after the SN explosion. In this case, the lines and the SN bump can not be seen in the same 
events. 

In our model, the disk was in place before the GRB occurred, so it may produce a high 
level of pre- GRB activity of the source. 

Different values of the ionization parameter, £, could produce the different reflection 
spectra. When £ ~ 10 2 , the spectra will show luminous lines from light metals and a de- 
pressed K a iron line; while 10 3 < £ < 10 5 , a luminous iron line will be observed in the spec- 
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tra(Lazzati et al. 2002). In our model the fallback disk with different properties, such as 
different of number density of the electrons in the disk surface and the outermost radius, will 
have different ionization parameter. In this paper, £ ~ 2.5 x 10 3 , so luminous K a iron line can 
be observed in the spectra. 

Vietri et al.(1999) has suggested a thermal model that a relativistic fireball associated with the 
GRB might hit the pre-GRB supernova remnant within ~ 10 3 s and heat the ejecta to T~ 10 7 - 10 8 
K. At such temperature the plasma emission shows thermal bremsstrahlung emission as well as 
iron line emission. In their model the thermal bremsstrahlung and recombination continuum from 
the thermal disk can account for the re-burst observed in GRB970508 and GRB970828. While in 
our model, the delayed injection energy from the central engine after the main burst, more than the 
energy of the main burst accounts for the re-burst and produces the iron line emission. 

Our model is also different from the decaying magnetar model in which Rees & Meszaros(2000) 
suggested that iron line could be attributed to the interaction of a continuing but decaying post- 
burst relativistic outflow from the central engine with the progenitor stellar envelope at distances 
less than a light-hour. In their model bumps should be found in less than several hours after the 
GRBs(Gao & Wei 2004). 

In conclusion, We suggest that the delayed injection energy that causes the re-burst in the 
GRB afterglow, illuminates the fallback disk which is formed after the supernovae explosion, pho- 
tonionizes the fallback disk, then produces iron line feature. This scenario can well explain the 
production of the re-burst and the emission lines, that can be tested by the observations of SWIFT 
satellite in the near future. 
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reflected iron photon 




Fig. 1. — This is just a cartoon picture of the geometry of the fallback disk. The central engine is 
a magnetar or a black hole. 



